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Abstract 

The proton-dissociative diffractive photoproduction of J ftp mesons has been 
studied in ep collisions with the ZEUS detector at HERA using an integrated 
luminosity of 112 pb" 1 . The cross section is presented as a function of the photon- 
^ ■ proton centre-of-mass energy and of the squared four-momentum transfer at the 

proton vertex. The results are compared to perturbative QCD calculations. 
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1 Introduction 



Photoproduction of vector mesons (VMs) is usually thought of as a process where the 
photon fluctuates into a qq state, which then interacts with the proton and becomes a 
VM. If the spatial configuration of the qq state is large, its interaction with the proton 
is soft in nature and is usually described by Regge theory pQ together with the vector 
dominance model [2]. This applies to exclusive photoproduction of the light VMs p, u and 
(f) (see Ivanov, Nikolaev and Savin [3] for a recent review). For heavy VMs, the qq pair 
is squeezed into a small configuration and perturbative QCD (pQCD) [I] can be applied. 
In exclusive photoproduction of J/tjj, 7 p — > J/iJj p, the mass of the J/ip provides a hard 
scale at the photon vertex and the small-size qq pair interacts through a two-gluon ladder 
with partons in the proton. If the four-momentum-transfer squared at the proton vertex 
is small, |i| < 1 GeV 2 , and the proton stays intact, the cross section is predicted to fall 
exponentially with \t\. 

When \t\ increases, \t\ > 1 GeV 2 , the dominant process is that where the proton dissociates 
into a low-mass nucleon state Y, 

jp-^J/ip Y. (1) 

At large \t\ values, the cross section is expected to have a power-law decrease with \t\ [5HH]. 
In addition, J/ip photoproduction at large \t\ is a two-scale process in which the large 
mass of the heavy VM is the hard scale at the photon vertex and t is the hard scale at the 
proton vertex. At high photon-proton centre-of-mass energies, W, this process should be 
sensitive to BFKL [9] dynamics. 

This paper contains results for the kinematic range 30 < W < 160 GeV and 2 < \t\ < 
20 GeV 2 , which is larger than for the previous ZEUS measurement [10]. The sample under 
study also represents more than a five-fold increase in integrated luminosity. 

2 Experimental set-up 

This analysis is based on data collected with the ZEUS detector at HERA in 1996—2000. 
In those years HERA operated with an electron^ beam energy of 27.5 GeV and a proton 
beam energy, E p , of 820 GeV (1996-1997) and 920 GeV (1998-2000). The data sample 
corresponds to an integrated luminosity of 112 pb -1 , 36 pb -1 with E p = 820 GeV and 
76 pb^ 1 with E p = 920 GeV. 

A detailed description of the ZEUS detector can be found elsewhere fl lj. A brief outline 
of the components that are most relevant for this analysis is given below. 

1 Electrons and positrons are both referred to as electrons in this paper. 
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Charged particles were tracked in the central tracking detector (CTD) [12], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consisted 
of 72 cylindrical drift chamber layers, organised in 9 superlayers covering the polar-angle^ 
region 15° < 6 < 164°. The transverse-momentum resolution for full-length tracks was 
<j(Pt)/pt = 0.0058pT©0.0065©0.0014/p r , withp T in GeV. The high-resolution uranium- 
scintillator calorimeter (CAL) [13] consisted of three parts: the forward (FCAL), the 
barrel (BCAL) and the rear (RCAL) calorimeters. Each part was subdivided transversely 
into towers and longitudinally into one electromagnetic section (EMC) and either one (in 
RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). The smallest subdivision 
of the calorimeter is called a cell. The CAL energy resolutions, as measured under test- 
beam conditions, were a(E)/E = 0.18 /y/~E for electrons and a(E)/E = 0.35/v^ for 
hadrons, with E in GeV. 

The muon system [T4] consisted of tracking detectors (forward, barrel and rear muon 
chambers: FMUON, B/RMUON), which were placed inside and outside a magnetised 
iron yoke surrounding the CAL. The inner chambers, F/B/RMUI, covered the polar 
angles from 10° to 34°, from 34° to 135° and from 135° to 171°, respectively. 

The luminosity was determined from the rate of the bremsstrahlung process ep — > cyp, 
where the photon was measured by a lead-scintillator calorimeter [TS] located at Z — 
-107 m. 



3 Kinematics and reconstruction 

The proton-dissociative J/ip production process in ep interactions, 

e(k)p(P) ->■ e(k')J/i/j(v)Y(P'), 

is illustrated in Fig. [TJ The signature of these events consists of two oppositely charged 
muons from the J/ip decay and of the remnant of the dissociated proton. In the case of 
photoproduction, the beam electron is scattered at small angles and escapes undetected 
down the beampipe. 

The variables k, k', P, P' and v are the four-momenta of the incident electron, scattered 
electron, incident proton, diffractive nucleonic system Y and J/ifj, respectively. The four- 
momentum of the exchanged photon is denoted by q. The kinematic variables are the 
following: 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. 
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• Q 2 = —q 2 = —(k—k') 2 , the negative squared four-momentum of the exchanged photon; 

• W 2 = (q + P) 2 , the squared centre-of-mass energy of the photon-proton system; 

• t = (P — P') 2 = (q — v) 2 , the squared four-momentum transfer at the proton vertex; 

• y = (P ■ q)/(P ■ k), the fraction of the electron energy transferred to the photon in the 
rest frame of the proton; 

• z = (P ■ v)/(P ■ q), the event inelasticity, i.e. the fraction of the virtual photon energy 
transferred to the J /if) in the proton rest frame. 

The dissociated proton either escapes undetected down the beampipe or deposits only a 
part of its energy in the CAL and hence the mass of the proton remnant, My, cannot 
be measured precisely. However, My is related to other kinematic variables through 
M 2 = W 2 (l -z)- \t\. 

The following angles are used to describe the decay of the J /if) (see Fig. [2]): 

• $, the angle between the electron-scattering plane and the vector-meson plane, in the 
photon-proton centre-of-mass frame; 

• 9h and (f>h, the polar and azimuthal angles of the positively-charged decay particle in 
the helicity frame. Here, the helicity frame is the J /if) rest frame and the quantisation 
axis is the meson direction in the photon-proton centre-of-mass system. The polar 
angle, 9h, is defined as the angle between the direction of the positively charged decay 
particle and the quantisation axis. The azimuthal angle, <frh, is the angle between the 
decay plane and the vector-meson production plane. 

In this study, photoproduction is characterised by the non-observation of the scattered 
electron. Thus, Q 2 ranges from the kinematic minimum, Q^ in = m 2 y 2 /(l — y) ps 
l(r 7 GeV 2 , where m e is the electron mass, up to Q^ax ~ 1 GeV 2 , the value at which 
the scattered electron becomes observable in the CAL. Since the mean Q 2 is small, 
{Q 2 ) ~ 5 ■ 10~ 5 GeV 2 , it was neglected in the reconstruction of the other kinematic vari- 
ables. 

The variable t can be expressed as t ~ — p\, where pt is the transverse momentum of 
the produced vector meson in the laboratory frame. The variable W is calculated as 
W 2 ~ 2E P (E — pz)j/ip, where E is the energy and pz is the longitudinal momentum of 
the vector meson. The quantities (E — pz)j/ip and t were reconstructed using only the 
measured momenta of the VM muon decay particles. 

The inelasticity z was computed from z — (E — Pz)j/4>/ XX ~ Pz), where ^2(E — pz) = 
(E — pz)j/4> + ^2(E — pz)had and ^2(E — pz)h&d is reconstructed by summing over all the 
CAL energy deposits (larger than 300 MeV) not associated with the J/tp candidate. 



3 



4 Event selection 



The events were selected online by the ZEUS three- level trigger system [TTJ[T6]. The 
events were required to have at least one track in the CTD. At least one track had to 
point towards a CAL energy deposit compatible with a minimum ionising particle as well 
as a signal in the inner muon chambers. 

The following was required offline: 

• no scattered electron observed; 

• two tracks with opposite charge pointing to a primary vertex with | Z vertex | < 50 cm; 

• both tracks well reconstructed, i.e. traversing at least three superlayers in the CTD, 
including the innermost layer; 

• each track associated with a distinct CAL energy deposit within a radius of 30 cm; 

• azimuthal angle between the two tracks associated with the two muon candidates less 
than 174° in order to reject cosmic-ray events; 

• invariant mass of the two tracks, which were assigned a /i mass, in the range 2.6 < 
M MM < 3.5 GeV. 

Events were required to be in a kinematic range where the properties of the final state 
particles were properly measured and the acceptance was well defined. This was satisfied 
for 2 < \t\ < 20 GeV 2 and 30 < W < 160 GeV. The cut of \t\ > 2 GeV 2 also significantly 
reduced the background from the exclusive process. A cut of z > 0.95 was applied to 
suppress non- diffract ive background. This cut also restricted the invariant mass of the Y 
system to M Y < 30 GeV. 

The energy range 30 < W < 40 GeV was mainly populated by events triggered by the 
FMUON detector, while the range 40 < W < 160 GeV was dominated by B/RMUON- 
triggered events. The FMUON-triggered sample was limited to the data collected in 
1996—1997 and covered the \t\ region up to 10 GeV 2 . 

After this selection procedure the number of observed di-muon events was 2817. 

5 Theoretical predictions 

The reaction 7 p — > J/ijj Y can be viewed as a three-step process. The photon fluctuates 
into a qq pair that scatters off a single parton in the proton by the exchange of a colour 
singlet. The scattered qq pair becomes a J/ip and the struck parton and the proton 
remnant together fragment into the system Y. In lowest-order QCD the colour singlet 
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exchanges a pair of gluons. In the leading logarithmic (LL) approximation, the process is 
described by the effective exchange of a gluon ladder. 

As stated in the introduction, the process under study has two scales. At the photon 
vertex, where the photon fluctuates into a qq pair, the size is fixed and determined by the 
J/ip mass. The second scale, controls the size of the system which emits the gluon 
ladder. 

In the region where the scale \t\ is smaller than Mj,^ (2 < \t\ < lOGeV 2 ), the momenta 
on the gluon ladder are still expected to be ordered and thus a DGLAP [17] approach 
is appropriate. A calculation in this kinematic region has been carried out by Gotsman, 
Levin, Maor and Naftali (GLMN) [18], using their screening correction formalism and 
evolving the gluon in a LL DGLAP mechanism. 

As \t\ increases, the BFKL mechanism is expected to dominate. The first LL BFKL 
calculations [SHZ] were made using the Mueller- Tang (MT) approximation [19], which is 
only good for very large rapidity intervals. Enberg, Motyka and Poludniowski (EMP) [5] 
do not use the MT approximation and provide a complete analytical solution in LL for the 
case of heavy quarks (the case for any quark mass is discussed elsewhere [2H1EI])- They 
use two different values of as as the pre-factor of the cross section and as the coupling 
relevant for the BFKL ladder. Enberg at al. [8] also present results of a non-leading 
(nonL) BFKL calculation. 

In addition, a recent QCD calculation by Frankfurt, Strikman and Zhalov (FSZ) [22|[23] 
is motivated by the QCD factorisation theorem for large \t\ rapidity-gap processes and by 
the correspondence to exclusive J/ip production at \t\ ~ 1 GeV 2 . In this QCD calculation 
in the triple-Pomeron limit, the W dependence of the cross section mainly depends on 
the gluon distribution of the proton. 

In all models, a non-relativistic approximation of the J/ip wave-function assuming equal 
sharing of longitudinal momenta between the quark and the anti-quark was used. The 
J/ip retains the helicity of the photon which means that s-channel helicity is conserved 
(SCHC). 

The DGLAP-motivated calculation predicts a mild W dependence of the cross section 
in the region of small \t\. In the region of larger \t\ the hard scale is chosen such that 
saturation is reached and thus the cross section is independent of W. The BFKL LL 
calculations predict a fast rise of the cross section with W which hardly depends on \t\. 
This is a unique feature of BFKL dynamics. The nonL BFKL model behaves in a similar 
way. In case of the FSZ parameterisation, the main energy dependence is provided by the 
behaviour of the gluon distribution. 

All calculations predict an approximate power-law t-dependence of the cross section of 
the form da/dt ~ \t\~ n , where the value of n may depend on the |t| range. 
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6 Monte Carlo and background evaluation 



The acceptance and the effects of the detector response were determined using Monte 
Carlo (MC) events. All generated events were passed through the standard ZEUS detector 
simulation, based on Geant 3.13 |24j . the ZEUS trigger-simulation package and the same 
reconstruction and analysis programs as used for the data. 

6.1 The process ep — >• e J/i(y Y 

The process ep —> e J ftp Y was modelled using the Epsoft generator J22H2B]. The 7p 
interactions were simulated assuming the exchange of a colourless object which couples 
to the whole proton, which subsequently fragments into a state Y. The particle multi- 
plicities and the transverse momenta of the hadrons in the final state Y were simulated 
using parameterisations of pp data, while the longitudinal momenta were generated with 
a uniform rapidity distribution. The differential cross-section da/dt was reweighted to 
obtain the shape observed in data. The assumption of s-channel helicity conservation was 
applied. 

The differential cross section in M y have the form da/dMy oc (My)"^^'. The measured 
z distributions in \t\ and W bins were used to determine the \t\ and W dependence 
of the function /3(t,W). For each bin, a single value of the function was extracted 
by using a x 2 minimisation method. The results were parameterised in the form of 
P(t,W) = (ty/iy ) a52±an exp((0.08±0.07) + (-0.14±0.03)|t|), with W = 95 GeV. This 
parameterisation was used in all further studies. 

6.2 Evaluation of background 

The main sources of background were the non-resonant QED 77 processes, misidenti- 
fied pion production and resonant background produced through the decay of the ip(2S) 
meson. 

The non-resonant background due to the QED Bethe-Heitler di-muon production, ep — > 
e/j, + /j,~Y, was simulated using the Grape-Dilepton 1.1 generator [27] . The background 
from 77 — > fi + fi~ events was estimated in each bin by normalising to the luminosity of 
the data. The contribution of this background increased with \t\ from 6 to 10%. 

The ■0(2iS f ) background was estimated using the Dipsi generator [2B]- This background 
was dominated by the processes tf> — > J/ip7i°7i° (BR=(17.51 ± 0.34)%) and ip — > /i + /i~ 
(BR=(0.76±0.08)%). It amounted to about 1% and 0.1%, respectively. The contribution 
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from ip — > J/ip7r + 7T (BR=(33.1 ± 0.5)%) was strongly suppressed by selection criteria, 
i.e. a requirement of two-track events. 

The background from exclusive J/ip production was found [10] to be 5% for 2 < \t\ < 
3GeV 2 . For |£| > 3GeV 2 , it was found to be consistent with zero. All background 
processes were subtracted bin-by-bin. 

7 Systematic uncertainties 

The systematic uncertainties were determined by varying the selection cuts and modi- 
fying the analysis procedure. Their effects on the integrated cross section are given in 
parentheses: 

• the cut on Z vertex was changed by ±10 cm (+1.8%, —0.5%); 

• the mass window was changed to 2.8 — 3.4 GeV (+0.2%); 

• instead of using the MC to subtract the background bin-by-bin, it was fitted with a 
polynomial function and statistically subtracted (+1.5%); 

• the minimum energy of the CAL energy deposit included for the evaluation of the z 
variable was varied by ±100 MeV (+0.3%, -1.4%); 

• the strategy of matching energy deposits to the decay tracks was changed. Instead 
of matching every object within 30 cm from the track, only one island within this 
distance was matched to the track (—0.7%); 

• the uncertainty of the muon acceptance, including the detector, the trigger and the 
reconstruction efficiency, was obtained from a study [29] based on an independent 
dimuon sample (±6.3%); 

• the uncertainty on the acceptance due to modelling of the hadronic final state in the 
EPSOFT MC was estimated by varying the parameter (3 within its errors (±2%). 

The overall systematic uncertainty was determined by adding all the individual uncer- 
tainties in quadrature. The uncertainty on the luminosity measurement, 2%, was not 
included. 

8 Results 

8.1 The J/?/> signal 

The invariant-mass distribution of the fi + fi~ pairs is presented in Fig. [31 A clear peak at 
the J /if) mass is observed with very little non-resonant background. 
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The distributions of the kinematic variables \t\, W, z, <ph and cos8h are shown in Fig. HJ 
The MC distributions of the J/ip events are shown as well as the QED background. The 
overall agreement between data and MC is good. 



8.2 Determination of photon-proton cross-section 

The ep cross section was determined by subtracting the background from the data, correct- 
ing for the acceptance, using the branching ratio for the muon channel decay (5.88±0.10%) 
and using the measured luminosity. 

Photon-proton cross sections were extracted from the ep cross sections by using photon 
flux factors. The flux factors [30] generated at the leptonic vertex relate the ep and the 
7P cross sections by 

where IV is the effective photon flux. The cross sections for different beam energies were 
averaged using the corresponding luminosities. 



8.3 \t\ dependence 

Differential ^p cross section da/dt for proton-dissociative J/ip photoproduction was mea- 
sured in the kinematic region 30 < W < 160 GeV, 2 < \t\ < 20 GeV 2 and z > 0.95. 

The differential cross section as a function of \t\ is shown in Fig. [5] and listed in Table [TJ 
The cross section falls steeply with \t\. The data cannot be described in the whole \t\ 
region by one exponential function of the form ~ e -6 '*', where b is a constant. Neither 
does a single power-law dependence of the form \t\~ n , where n is a constant, fit the data. A 
good fit can, however, be obtained by fitting two \t\ ranges separately, giving n = 1.9 ±0.1 
for 2 < \t\ < 5 GeV 2 and n = 3.0 ± 0.1 for 5 < \t\ < 20 GeV 2 . Note that a good fit can 
also be obtained to a quadratic exponential function g-^l+cl*! 2 . 

The differential cross-section da/dt as a function of \t\ is shown again in Fig. [6], together 
with the HI data [31], and compared with different theoretical models. The GLMN LL 
model gives a good description of the data up to about \t\ — 5 GeV 2 , but falls off slower 
than the data up to the region where the calculation is valid (\t\ < 10 GeV 2 ). The EMP LL 
prediction, using as = 0.205 in the pre-factor and as = 0.16 in the BFKL evolution, lies 
below the data in the whole range of \t\. The FSZ results are shown for a calculation using 
a Pomeron trajectory with intercept of 1.1 and a slope of 0.005 GeV~ 2 . Similar results 
are obtained with a Pomeron intercept of 1.0. The CTEQ6M parameterisation [32] of the 
parton density functions is used. The FSZ calculation describes the data well up to \t\ of 
about 12 GeV 2 but falls-off too steeply at larger |t| values. 
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8.4 W dependence 



In the Regge formalism, the differential cross section can be expressed as 

da/dt = F{t)W AMt) - l \ (2) 

where F(t) is a function oft and a P (t) is the effective Pomeron trajectory. This expression 
is usually used for exclusive reactions, but has been used also for the case where M Y is 
integrated over [TQlEl] . By studying the W dependence of da/dt at fixed t, the Pomeron 
trajectory can be determined. 

The W dependence of the differential cross section da/dt for eight fixed t values is shown 
in Fig. [7] and listed in Table [21 At each t value the cross section is parameterised as 
a ~ W 5 and the lines in the figure are the result of these fits. The values of a P can be 
obtained at each t value through 

a P = (5 + 4)/A (3) 

and are shown in Fig. [HJ The values of 6 and of a P are listed in Table [3j A linear fit of 
the form 

a P (t) = a P (0) + a' p ■ t (4) 

yields an intercept 

a P (0) = 1.084 ± 0.031 (stat.)l^(syst.), (5) 

and a slope 

a' P = -0.014 ± 0.007(stat.)i^(syst.). (6) 

The value of the intercept is consistent with that of the so-called "soft" Pomeron [33] 
(1.0808). The slope is different from that of the "soft" Pomeron [34] (0.25 GeV~ 2 ), but is 
consistent with the predictions of the BFKL Pomeron [351136] . 

The 7p cross section as a function of W was measured in four bins of|£|:2<|£|<3 GeV 2 ; 
3 < \t\ < 5 GeV 2 ; 5 < \t\ < 10 GeV 2 and 10 < \t\ < 20 GeV 2 for 30 < W < 160 GeV. The 
cross-section values are shown in Fig. [9] and summarised in Table HI The HI data [31] for 
the \t\ bin of 5 to 10 GeV 2 are also shown. A clear rise with W is seen in all the four \t\ 
regions. Also shown in the figure are the predictions of the models used in the comparison 
with da/dt. The DGLAP-based GLMN LL calculation agrees well with the data in the 
first two \t\ bins, but fails to describe the rise with W for \t\ > 5 GeV 2 . The other two 
calculations, EMP LL and FSZ, predict a W dependence which is too steep in all the \t\ 
ranges presented in the analysis. 
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8.5 Decay angular distributions 



The angular distributions of the J/ip decay provide information about the photon and 
J ftp polarisation states. The normalised two-dimensional angular distributions can be 
written in terms of spin density matrix elements, r, as: 

^l + rE)-±(3r£-l)cos 2 fc + 



ad cos B h d(f) h 4tt \2 y uuy 2 

+v / 2Re{r° 4 } sin 26 h cos <p h + r°i 1 sin 2 6 h cos 2<f) h ) . (7) 

The one- dimensional distributions result from the integration over 9^ or <ph and are ex- 
pressed as : 

ocl + r ° 4 + (l-3r ° 4 )cos 2 ^ (8) 



d cos Oh 
and 

^cxl + r? 4 . 1 cos20 h . (9) 
d(p 

The spin density matrix element Tq 4 represents the probability that the produced J/ip has 
helicity zero, Re{r° 4 } is proportional to the single-flip amplitude and r^_ 1 is related to the 
interference between non-flip and double-flip amplitudes. If SCHC holds, the J/ip retains 
the helicity of the almost real photon and all the three matrix elements are expected to 
be zero. 

The distributions of cos 9h and <J>h after background subtraction and acceptance corrections 
in four \t\ bins and for 30 < W < 160 GeV are shown in Fig. [TD1 They were fitted using 
formulae ([ED and (Q. The r^ 4 , Re{r| 4 } and rjij spin density matrix elements were 
extracted from a two-dimensional x 2 minimisation fit using Eq. (j7j) and are summarised 
in Table [5] and shown in Fig. [TTJ The HI data [31] are also shown. The results for r§ 4 
and r^_ x are compatible with zero. The values obtained for Rejr'j' 4 } are not compatible 
with zero for \t\ < 10 GeV 2 , contrary to the expectation from SCHC. 

The measurements of the present analysis, the differential cross section as a function of 
\t\, the W dependence of the cross section and the density matrix elements, are in good 
agreement with those of the HI collaboration [3T] within the common kinematic region. 



9 Summary 

Proton-dissociative J / ip production was measured at HERA in the photoproduction regime 
in the kinematic region 30 < W < 160 GeV, z > 0.95 and 2 < \t\ < 20 GeV 2 . 
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The \t\ dependence of the differential cross section, da/d\t\, is found to be approximately 
power-like, ~ \t\~ n , with the power n increasing with \t\. 

The effective Pomeron trajectory was derived from a measurement of the W dependence 
of the cross section at fixed t values. The value of the slope of the trajectory is compatible 
with zero. It is consistent with the predictions of the BFKL Pomeron but different from 
the slope of the "soft" Pomeron. 

The cross-section a(jp —> J/ip Y) rises significantly with W in each \t\ bin. The t and W 
dependence of the cross section were compared to several theoretical calculations. The 
DGLAP-motivated GLMN LL [18] calculation can describe the behaviour of the data, 
both in t and in W, up to \t\ = 5 GeV 2 . The BFKL-motivated EMP LL [S] calculation fails 
to describe the data in the kinematic region of the present measurement. The FSZ [22|l23] 
calculation describes the t dependence of the cross section only up to \t\ = 12 GeV 2 and 
fails to reproduce the W dependence. 

The spin density matrix elements of the J/i/j, and are consistent with zero, as 
expected from s-channel helicity conservation. The values obtained for Re{r°Q} are not 
compatible with zero for \t\ < 10 GeV 2 , contrary to the expectation from SCHC. 
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|t| bin 
( GeV 2 ) 


(1*1) 
( GeV 2 ) 


da/d\t\ 
(nb/GeV 2 ) 


2.0 3.0 
3.0 - 4.0 


2.5 
3.5 


6.05 ± 0.23 toll 
3.14 ± 0.18 toil 


4.0 - 5.0 


4.4 


1.91 ± 0.13 toil 


5.0 - 6.5 


5.6 


0.97 ± 0.08 tlH 


6.5 - 8.0 


7.2 


0.41 ± 0.05 t°o° i 


8.0 - 11.0 


9.3 


0.21 ± 0.02 to° i 


11.0 - 14.0 
14.0 - 20.0 


12.4 
16.5 


0.07 ± 0.01 t^j 
0.05 ± 0.01 t°o° \ 



Table 1: Differential cross-section do~/d\t\ as a function of \t\ for 30 < W < 
160 GeV and z > 0.95. The first uncertainty is statistical and the second is sys- 
tematic. 
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\t\ bin 
( GeV 2 ) 



2.0-2.5 



2.5-3.0 



3.0-4.0 



( GeV 2 ) 



2.2 



2.7 



3.4 



W bin 

(GeV) 



30 
50 
70 
80 
90 
100 
110 
120 
130 



50 
70 
80 
90 
100 
110 
120 
130 
160 



30 
50 
70 
80 
90 
100 
110 
120 
130 



50 
70 
80 
90 
100 
110 
120 
130 
160 



30 
50 
70 
80 
90 
110 
130 



50 
70 
80 
90 
110 
130 
160 



(W) 
(GeV) 



39.5 
59.5 
75.0 
84.8 
95.0 
105.1 
114.9 
125.1 
144.0 



39.6 
59.3 
74.8 
85.0 
95.0 
104.8 
114.7 
124.9 
144.5 



39.7 
59.5 
74.8 
84.8 
99.8 
119.6 
143.5 



da/dt 
nb/GeV 2 ) 



3.82 ± 0.59 toil 
6.65 ± 0.61 tof 5 
8.52 ± 1.10 t° f 8 
5.95 ± 0.91 tofi 



7.10 ± 1.08 



h0.57 
-1.18 



7.51 ± 1.13 t f 



+0.99 
-0.58 
+0.98 
-0.78 

6.74 ± 0.93 tlf 3 



6.83 ± 1.20 
8.79 ± 1.52 



2.47 ± 0.52 



+0.60 
-0.31 



3.93 ± 0.50 tof 3 

+0.74 
-0.91 



5.43 ± 0.91 



4.52 ± 0.84 tojl 
4.89 ± 0.93 toil 
4.81 ± 0.95 t°of 3 
6.02 ± 1.26 toll 
6.34 ± 1.33 tUt 
7.00 ± 1.17 t° f 8 



2.06 ± 0.34 toll 
2.65 ± 0.30 toll 
3.33 ± 0.53 toil 
3.17 ± 0.53 toll 
3.38 ± 0.45 toil 
3.72 ± 0.54 toil 

+0.49 
-0.38 



3.26 ± 0.52 



Table 2: Differential cross section as a function 
z > 0.95. The first uncertainty is statistical and the 



of W in eight t bins and for 
second is systematic. 
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|t| bin 


-t 


W bin 


(W) 


da/dt 


( GeV 2 ) 


( GeV 2 ) 


(GeV) 


(GeV) 


( nb/ GeV 2 ) 






30 - 50 


39.5 


0.95 ± 0.23 toil 






50 - 70 


59.8 


1.51 ± 0.22 to'17 






70 - 80 


74.8 


1.97 ± 0.38 + no9 


4.0-5.0 


4.5 


80 - 90 


85.1 


1.98 ± 0.39 






90 - 110 


99.6 


1.71 ± 0.29 +0 n it 

— U.oU 






110 - 130 


119.7 


2.44 ± 0.44 +0 n ii 

w^" —0.29 






130 - 160 


144.5 


2.82 ± 0.51 toil 






30 - 70 


48.1 


0.45 ± 0.15 tom 






50 - 70 


59.6 


0.79 ± 0.14 torn 


5.0-6.5 


5.7 


70 - 90 


79.9 


1.00 ± 0.15 toll 






90 - 110 


100.0 


1.07 ± 0.19 +0 ]l 






110 - 160 


134.1 


1.20 ± 0.17 to.ll 






30 - 70 


47.7 


0.17 ± 0.11 






50 70 


60.0 


0.28 ± 0.08 tom 


6.5-8.0 


7.2 


70 - 90 


79.4 


0.22 ± 0.07 tom 






90 - 110 


99.8 


0.43 ±0.11 + Hfi 






110 - 160 


133.6 


0.64 ± 0.13 t'oll 






50 - 70 


59.5 


0.12 ± 0.04 to l 


8.0-11.0 


9.2 


70 90 
90 - 110 


80.4 
99.6 


0.24 ± 0.05 tof 3 
0.21 ± 0.06 torn 






110 - 160 


133.4 


0.26 ± 0.05 t°o°oi 






50 - 70 


59.8 


0.04 ± 0.01 tom 


11.0-20.0 


14.2 


70 - 90 
90 - 110 


79.3 
100.0 


0.05 ± 0.01 tofi 
0.07 ± 0.02 +°;° 2 






110 - 160 


133.9 


0.07 ± 0.02 to°ol 



Table 2 (Continuation): Differential cross section as a function of W in eight t bins 
and for z > 0.95. The first uncertainty is statistical and the second is systematic. 
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|t| bin 
( GeV 2 ) 


-t 

( GeV 2 ) 


5 




2.0 


- 2.5 


2.2 


0.38 


± 


0.10 


+0.11 
-0.06 


1 


10 


± 





03 


+0.03 
-0.02 


2.5 


- 3.0 


2.7 


0.70 


± 


0.15 


+0.10 
-0.10 


1 


18 


± 





04 


+0.02 
-0.02 


3.0 


- 4.0 


3.4 


0.39 


± 


0.13 


+0.08 
-0.07 


1 


10 


± 





03 


+0.02 
-0.02 


4.0 


- 5.0 


4.5 


0.72 


± 


0.18 


+0.08 
-0.07 


1 


18 


± 





05 


+0.02 
-0.02 


5.0 


- 6.5 


5.7 


0.70 


± 


0.21 


+0.12 
-0.12 


1 


18 


± 





05 


+0.03 
-0.03 


6.5 


- 8.0 


7.2 


1.38 


± 


0.46 


+0.37 
-0.24 


1 


35 


± 





12 


+0.09 
-0.06 


8.0 


- 11.0 


9.2 


0.78 


± 


0.38 


+0.06 
-0.12 


1 


20 


± 





09 


+0.02 
-0.03 


11.0 


- 20.0 


14.2 


0.82 


± 


0.44 


+0.34 
-0.08 


1 


21 


± 





11 


+0.09 
-0.02 



Table 3: The values of the parameters 5 and the effective Pomeron trajectory 
ctp for eight fixed t values. The first uncertainty is statistical and the second is 
systematic. 
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\t\ bin 


-t 


W bin 


(W) 




da/dt 


( GeV 2 ) 


( GeV 2 ) 


( GeV) 


(GeV) 


(nb/GeV 2 ) 






ou 




^n 
ou 


39.5 


'i 9^ 
o.zo 


4- 


n An +0.43 
U - 4U -0.33 






ou 




7n 
( u 


59.4 


o.oo 


_1_ 


n 4n +0.31 
U - 4U -0.34 






( u 




sn 
oU 


74.9 


( .16 


i 

it 


n 79 +0.74 
u - ' z -0.49 






sn 
ou 




on 
yu 


84.9 


K 'in 
o.ou 


_1_ 


n fi9 +0.47 
U -° 2 -0.50 


2-3 


2 5 


90 




100 


95.0 


6 04 

W . Wt: 


-£ 


71 

u - ' 1 -0.45 






1 nn 

1UU 




1 1 n 

1 1U 


105.0 


fi 98 
u.zo 


4- 


n 7 cr +0.57 
u - ' -0.37 






1 1 n 

J. J.U 




1 2n 


114.8 




4- 


86 + - 57 
u.so _ 41 






1 2D 




1 ^n 


125.1 


7.52 


± 


1 oo +0 - 77 

i - uu -0.60 






130 


— 


160 


144.3 


6.80 


± 


72 +0 - 71 

u - ' z -0.57 






ou 




^n 
ou 


39.7 


3.10 


± 


42 +0 - 21 

u -^ z -0.20 






^n 

ou 




7n 

( u 


59.6 


4.23 


± 


n vj +0.22 

u -°' -0.23 






7n 

i u 




sn 

ou 


74.9 


5.28 


± 


65 +0 - 31 

U - Ud -0.36 






sn 

OU 




on 
yu 


84.9 


5.22 


± 


67 +0 - 48 

u -°' -0.45 


3-5 


3.8 


90 




100 

X WW 


94.9 


4.29 


± 


63 +0 - 41 

U.UO _ 26 






100 




110 


105.0 


6.42 


± 


n Q-i + - 68 
u.jo _ 42 






110 




120 


115.0 


5.79 


± 


92 +0 - 39 

U - JZ -0.32 






120 




130 


125.0 


6.55 


± 


1 07 +0 - 54 

-0.53 






130 




160 


143.9 


6.40 


± 


75 +0 - 77 

u - '° -0.63 






30 




70 


53.9 


1.69 ± 0.20 toil 


5-10 


6.7 


70 
90 




90 
110 


79.9 
99.9 


2.36 ± 
2.69 ± 


0.28 „°; 2 ° 
0.35 ±8^| 






110 




160 


133.8 


3.61 ± 0.36 ±g;|| 






50 




80 


63.9 


0.50 ± 0.108 to° 7 


10-20 


13.3 


80 




120 


98.8 


0.72 ± 0.132 t° f 2 






120 




160 


139.8 


0.84 ± 0.188 t°oi 4 



Table 4: 27ie cross-section <r(7p — > J/?/> V) as a function of W in four \t\ bins 
and for z > 0.95. The first uncertainty is statistical and the second is systematic. 
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\t\ bin 
(GeV 2 ) 


<l*l> 
( GeV 2 ) 


r 04 
' 1-1 


r 04 
'00 


£>e{r? 4 } 


2 3 
3-5 
5-10 
10 - 20 


2.5 
3.8 
6.7 
13.3 


0.005 ± 0.064 ±g;g|| 
— 0.206 ± 0.072 

0.003 ± 0.106 ±g;gfg 
— 0.164 ± 0.240 iftffjj 


0.090 ± 0.088 ±g;g$? 
— 0.030 ± 0.100 t otl 
-0.033 ± 0.147 IK 
-0.259 ± 0.328 


0.117 ± 0.061 ±g;gfg 
0.197 ± 0.068 ±S:os5 
0.154 ± 0.088 ±g;g^f 
— 0.153 ± 0.172 tow 



Table 5: The spin density matrix elements for 30 < W < 160 GeV and z > 0.95. 
The first uncertainty is statistical and the second is systematic. 
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Y(P') 



Figure 1: Schematic diagram of proton-dissociative J/ip production in ep inter- 
actions, ep — > eJ/ip Y . 

yp centre-of-mass 




electron-scattering plane production plane 



JA|/ direction 
in the yp 
centre-of-mass system 




J/\|/ rest frame 

Figure 2: Angles used to analyse the helicity states of the J/ifi, see text. 
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Figure 3: T/ie invariant-mass spectrum for pairs in the range 30 < 

W < 160 GeV , 2 < \t\ < 20 GeV 2 and z > 0.95. Error bars represent only 
statistical uncertainties. The data are compared to the MC distributions. The 
hatched histogram represents the ep — > efi + fi~Y background as simulated by the 
Grape MC. The solid-line histogram represents the sum of J/ip and background 
MC events. 
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Figure 4: Comparison between the data and MC distributions in the range 
30 < W < 160 GeV, 2 < \t\ < 20 GeV 2 and z > 0.95 for a) \t\, b) W, c) z, d) <j) h , e) 
cos Oh- Error bars represent only statistical uncertainties. The hatched histograms 
represent the ep — > efi + fi~Y background as simulated by the Grape MC. The 
solid-line histogram represents the sum of J/ip and background MC events. 
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ZEUS 




Figure 5: The \t\ dependence of the differential cross-section da/d\t\ for the 
process ^p — >■ J/ip Y at (W) = 81 GeV and z > 0.95. The inner bars correspond 
to the statistical uncertainties and the outer to the statistical and systematic un- 
certainties added in quadrature. The solid lines are the results of power fits to the 
form da/dt ~ \t\ n . 



23 



ZEUS 



CM 
> 

o 
O 

JO 



10 



t 

T3 



10 



• ZEUS 112 pb" 1 
□ H1 78 pb" 1 
EMP LL 



GLMN LL 



10 




<W> = 81 GeV 
z > 0.95 



2 4 6 8 10 12 14 16 18 20 



(GeV 2 ) 



Figure 6: The \t\ dependence of the differential cross-section da/d\t\ for the 
process ^p -> J/ip Y at (W) = 81 GeV and z > 0.95. The HI data, 50 < 
W < 150 GeV , J31\j are also shown. The inner bars correspond to the statistical 
uncertainties and the outer to the statistical and systematic uncertainties added in 
quadrature. The lines show the predictions of several calculations, referred to in 
the text. 
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Figure 7: The W dependence of the differential cross-section da/dt for the 
process ^p — > J/ip Y (z >0.95) at fixed \t\ values, as indicated in the figure. The 
inner bars correspond to the statistical uncertainties and the outer to the statistical 
and systematic uncertainties added in quadrature. The solid lines are the results of 
fits to the form da/dt ~ W s . 
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Figure 8: The effective Pomeron trajectory as a function of t. The inner 
bars correspond to the statistical uncertainties and the outer to the statistical and 
systematic uncertainties added in quadrature. The solid line is a fit of the form 
a F {t) = ajp(0) + a' p ■ t. The dashed line is an extrapolation to a P (0). 
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Figure 9: The W dependence for the process ^p — > J/ip Y (z >0.95) in jour dif- 
ferent \t\ bins. The HI data I31\j for the \t\ bin of 5 to 10 GeV 2 are also shown. The 
inner bars correspond to the statistical uncertainties and the outer to the statistical 
and systematic uncertainties added in quadrature. The lines show the predictions 
of several calculations referred to in the text. 
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Figure 10: TTie normalized distributions of <f> and cos(#) /or 30 < < 160 GeV 
and 2 > 0.95 in four bins of \t\. Error bars represent only statistical uncertainties. 
The lines represent the results of the fits according to formulae |2P and in the 
text. 
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Figure 11: Helicity spin density matrix elements a) b) Tqq and c) Rejr^Q} 
as a function of \t\ in the range 30 < W < 160 GeV and z >0.95. The HI 
data, 50 < W < 150 GeV , ISDj are also shown. The inner bars correspond to the 
statistical uncertainties and the outer to the statistical and systematic uncertainties 
added in quadrature. The solid lines show the expectation from SGHC. 
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